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Abstract. We present Ulysses and NEAR data from the detection of the short or intermediate duration (2 s) gamma-ray burst 
GRB 00030 1C (2000 March 1.41 UT). The gamma-ray burst (GRB) was localised by the Inter Planetary Network (IPN) and 
RXTE to an area of ~50 arcmin 2 . A fading optical counterpart was subsequently discovered with the Nordic Optical Telescope 
(NOT) about 42 h after the burst. The GRB lies at the border between the long-soft and the short-hard classes of GRBs. If 
GRB 00030 1C belongs to the latter class, this would be the first detection of an afterglow to a short-hard burst. We present 
UBRI photometry from the time of the discovery until 11 days after the burst. We also present JHK photometry obtained 
with UKIRT on 2000 March 4.5 UT (3.1 days after the burst). Finally, we present spectroscopic observations of the optical 
afterglow obtained with the ESO VLT Antu telescope 4 and 5 days after the burst. The optical light curve is consistent with 
being achromatic from 2 to 1 1 days after the burst and exhibits a break. A broken power-law fit yields a shallow pre-break 
decay power-law slope of ai = — 0.72±0.06, a break time of tbreak = 4.39±0.26 days after the burst, and a post-break slope 
of «2 = — 2.29±0.17. These properties of the light curve are best explained by a sideways expanding jet in an ambient medium 
of constant mean density. In the optical spectrum we find absorption features that are consistent with Fe II, C IV, C II, Si II and 
Lya at a redshift of 2.0404±0.0008. We find evidence for a curved shape of the spectral energy distribution of the observed 
afterglow. It is best fitted with a power-law spectral distribution with index /3 ~ —0.7 reddened by an SMC-like extinction law 
with Ay- ~ 0. 1 mag. Based on the Lya absorption line we estimate the H I column density to be log(N(H l))= 21 .2 ± 0.5. This 
is the first direct indication of a connection between GRB host galaxies and Damped Lya Absorbers. 



1. Introduction 



The discovery of the first X-ray afterglow (Costa et al. 1997) 
and optical counterpart (van Paradijs et al. 1997| ) to a long- 
duration gamma-ray burst (GRB) have led to a revolution in 
GRB research. The determination of a redshift of 0.835 for 



GRB 970508 (Metzger et al. 1997), and the subsequent de 



termination of redshifts of 13 bursts with a median redshift 
of ~1.0, have firmly established their cosm ologic al origin 
(Kulkarni et al. |2000aj This work; Bloom et al. |2000| ). 

The intriguing case of an association of the pecu- 
liar s uperno va SN1998bw with GRB 980425 (Galama 
et al. 1998) was the first indication of a possible con- 



nection with supernovae. Evidence for supernova signatures 
in the late lig ht cur ves of GRB 970228 (Reichart |1999 ; 
Galama et al. Il999|) and GRB 980326 (Castro-Tirado & 



Gorosabel 1999; Bloom et al. 1999 1 suggests that at least 
some long-duration GRBs may be related to the collapse of 
massive (> 25 M ) stars. Breaks in the power-law declines 
of GRB 99012 3 (Ku lkarni et al. |l999| ) and GRB 990510 
(Harrison et al. 1999 1 are interpreted as evid ence for colli- 
mated outflows ('jets') (see also Holland et al. 2000 ). Further 
evidence for th is collapsar + jet model (e.g., MacFadyen & 
Woosley |l999| ) comes from the light curve of GRB 980519 
which is best interpreted as a jet expa nding i nto a preexisting 
circumburst stellar wind (Jaunsen et al. 2001 ). 

The high-energy properties of GRBs show a bi-modal dis- 



tribution of burst durations (Kouveliotou et al. 1995) which, 



in the simplest scenario, may indicate the existence of binary 
compact mergers as the progenitors of the short-duration bursts 
(T90 < 2 s). Fro m an a nalysis of the Third BATSE Catalog, 
Mukherjee et al. ( 1998 ) have shown that, in addition to the 
short (Tgo < 2 s) and long (T90 > 5 s) classes, there may 
exist a third, intermediate soft-spectrum class of GRBs with 
duration 2 s < Tgo < 5 s. 

In this paper we report the discovery and subsequent 
observations and analysis of the afterglow of the short-to- 



intermediate duration GRB 00030 1C (Fynbo et al. 2000a) 



Sect. U reports the detection, IPN localisation and the high- 
energy data of the GRB obtained from Ulysses and NEAR. 
Sect. H describes the discovery of the optical counterpart and 
our subsequent optical and infrared observations. Sect. |] de- 
tails the optical and infrared photometry and Sect. || describes 
the VET spectroscopy. Sect. ^ describes the results obtained on 
the spectroscopy and spectral energy distribution and Sect. [7] is 
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devoted to the discussion and interpretation, with Sect. |8J pre- 
senting our conclusions. 

Throughout this paper, we adopt a Hubble constant of Ho = 
65 km s _1 Mpc -1 and assume f_ m = 0.3 and = 0.7. 



2. Detection and localisation of the gamma-ray 
burst 

GRB 00030 1C was recorded by the Ulysses GRB experiment 
and by the NEAR X-Ray/Gamma-Ray Spectrometer. Because 
this burst was relatively weak, it did not trigger the Ulysses 
Burst Mode, and the only data available from Ulysses is the 
Observation Mode 1 0.25 s resolution 25-150 keV light curve 
(Hurley et al. 1992b ). NEAR records the light curves of bursts 
in the 150-1000 keV energy range with 1 second resolution, 
but takes high-energy spectra only with 40 min resolution. 

Analysis of the Ulysses and NEAR relative tim- 
ing data yields an annulus centred at (a, <5)20oo = 
(20 h 34 m 7.56 s ,+20°32'19.62"), with a radius of 
57.520±0.083 degrees (at 3a full-width). This annulus 
intersected the error-box of the All-Sky Monitor (ASM) on the 
RXTE spacecraft, at near-right angles to create a composite 
localisation of a parallelogram of area 50 arcmin 2 (see Fig. ||). 

Since no high-energy spectra are available, we have esti- 
mated the peak fluxes and fluences for trial power-law spec- 
tra with indices between 1 and 4 using the Ulysses data. For 
a typical power-law index of 2, we find a 25-100 keV flu- 
ence of 2.1 x 10~ 6 ergcm~ 2 , and a peak flux over the same 
energy range, and over 0.25 s, of 6.3 x 10~ 7 erg cm -2 s _1 . 
The uncertainties in these numbers are partly due to the lack 
of a high-energy spectrum. For example, the fluence estimates 
range from 1.45 x 10~ 6 to 2.24 x 10 -6 erg cm -2 as the spec- 
tral index is varied from 4 to 1. The statistical uncertainty is ap- 
proximately 30%. From the NEAR data we estimate the 150- 
1000 keV fluence to be approximately 2 x 10~ 6 erg cm~ 2 . 

To date, the only GRBs with identified long-wavelength 
counterparts have been long-duration bursts. As measured by 
both Ulysses and NEAR, in the >25 keV energy range, the 
duration of this burst was approximately 2 s. (Note that the 
earlier estim ate of a 10 s duration of GRB 000301C by Smith 
et al. (2000) was based on the <10 keV energy range). Thus 
it falls in the short class of bursts, though it is consistent 
with belonging to the proposed intermediate class or the ex- 
treme short end of the distribution of long-duration GRBs 
(Hurley et al. |l992a| ; Mukherjee et al. |1998| ). Although we 
do not have any measurements of the high-energy spectra 
above 25 keV, it is possible to derive a crude estimate of the 
spectral index, and therefore the hardness ratio (the 100-300 
keV fluence divided by the 50-100 keV fluence), from the 
Ulysses and NEAR count rates. We obtain a hardness ratio 
of 2.7±0.6(cutoff)±30%(statisticalerror) from fitting a power- 
law, with the index as a free parameter, to the count rates from 
NEAR and Ulysses, assuming a range of cut-off energies. 

Fig. [j] shows the location of GRB 000301C in a hardness 
vs. duration plot. The contour plot contains 1959 GRBs for 
which data on fluence and duration were available in the Fourth 



BATSE GRB Catalog (revised) (Paciesas et al. 1999) and the 
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Fig. 1. A contour plot showing the duration-hardness 
(log(T 9 o)-log(HR 3 2)) distribution of BATSE bursts (fr om the 
Fourth BATSE GRB Catalog (revised) (Paciesas et al. |1999|) 



and the Current BATSE GRB Catalog). The triangle with an 
error-bar near the center of the plot represents GRB 000301C. 
Other symbols represent 10 other BATSE bursts with identi- 
fied counterparts for which data on fluence and duration are 
available. Triangles are bursts which have a break in their opti- 
cal light curves. Errors in the BATSE data are smaller than the 
symbol size. Contour levels scale linearly. The centroid in the 
lower left corner indicates the resolution. 

BATSE Current GRB Catalog]]. The symbols represent the 10 
GRBs included in the BATSE catalogs for which an afterglow 
has been identified, with GRB 00030 1C located near the center 
of the plot. Triangles are bursts where a break has been found 
in the optical light curve. From this sparse set of data there does 
not appear to be any marked difference in the distributions of 
bursts with, or without, an identified break. 

Of the 1959 BATSE bursts in Fig. |, the ratio between 
bursts with a duration of T90 > 2.0 s and with T90 < 2.0 s 
is 3:1. To date, at least 23 GRB optical afterglows have been 



discovered (Kulkarni et al. 2000a; Andersen et al. 2000; this 



work; Kloseet al. |2000[ Fynbo et al. |2001a ; Fynbo et al. [2001b ; 
Henden 2001; grb-webpage of J. GreinerQ). If the distribution 
of the 23 GRBs with identified counterparts follows the gen- 
eral BATSE distribution, one would expect that 17 ± 4 bursts 
were in the long class, and 6 ± 2 bursts were in the short class. 
However, GRB 00030 1C is the only GRB with a duration con- 
sistent with the short-duration class. The expected number of 
identified short burst counterparts is moderated by the strong 
selection bias caused by the technical difficulties of obtaining 
precise localisations for the short GRBs. 



3. Discovery and observations of the afterglow 



The IPN/RXTE error-box of GRB 00030 1C (Smith et al. 200C ) 



was observed with the 2.56-m Nordic Optical Telescope (NOT) 



Data on current GRBs are available through the BAT SE homepage 



http : / /www . batse . msf c . 
1 [http: 77 www . aip . de/ " 



. qov/batse/ 



jcg/ grbgen . html 




Fig. 2. A finding chart with the rectangular IPN-errorbox indi- 
cated along with the two initial ALFOSC pointings (ALFOSC 
FOV = 6.5 ')■ 



on 2000 March 3.14-3.28 UT (-1.8 days after the burst) 
using the Andalucfa Faint Object Spectrograph (ALFOSC). 
Comparing with red and blue Palomar Optical Sky Survey II 
exposures, a candidate Optical Transient (OT) was found at the 
positi on [a. (5)2ooo = (16 h 20 m 18.56 s , +29°26'36.1") (Fynbo 
et al. 2000a). A finding-chart of the IPN-errorbox and the two 
ALFOSC pointings used to cover the field can be seen in Fig. |[ 
The transient nature of the candidate was subsequently 
confirmed at optical and infrared wavelengths (Bernabei et 
al. [2000|; S tecklum et al |200C| ; Garnavich et al. |2000[ Veillet 
et al. |200C| ; Fynbo et al. [2000b| ; Kobayashi et al. ^00C| ). At the 
time of discovery the magnitude of the OT was R=20.09±0.04 
(see Sect. Q for a detailed discussion of the photometry and 
Fig. ^ for a finding chart and UBRI-images of the OT on 2000 
March 3 UT). 

We obtained subsequent optical observations using the 
NOT, the Antu telescope (UT1) of ESO's Very Large Telescope 
(VLT), the USNOFS 1.0-m telescope, the 2.2-m telescope at 
Calar Alto (CAHA) and the Wide Field Imager (WFI) on the 
ESO 2.2-m telescope. In addition we obtained near-infrared 
(NIR) data from the United Kingdom Infra-Red Telescope 
(UKIRT). 

The journal of our optical and NIR observations, including 
the derived magnitudes, is given in Table [IJ 

The OT was also observed from several other optical and 
infrared telescopes, and the counterpart was subsequently de- 



tected in the mm-band at 250 GHz (Bertoldi 200C ), and at ra- 
dio (8.46 GHz) wavelengths (Berger & Frail |200C[ )f|. Papers 
detailing the properties of the counterpart of GRB 00030 1C 



See the GCN Circulars Archive for further details 



(http: / /gcn.gsfc .nasa. go v / gen/ gen 3_ar chive .html) 
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Fig. 3. Upper panel: A contour plot showing a 3.5x3.5 arcmin 2 
region centred on the Optical Transient (OT) of GRB 00030 1C 
from the combined R-band image obtained at the NOT on 
2000, March 3 UT. The OT is located 6 arcsec east of the star 



marked A (following Garnavich et al. 2000 1. Lower panel: A 
region of 30x30 arcsec 2 centred on the OT from the combined 
U, B, R and I frames obtained on 2000, March 3.14-3.28 UT, 
about 42 hours after the burst. 



have been presented at radio and mm wavelengths (Berger et 
al. 



2000 1 and in the infrared (Rhoads & Fruchter 



2001), and 



optical bands (Masetti et al. 2000; Sagar et al. 200C). 



Table 1. Journal of our observations of GRB 000301C with 
NOT+ALFOSC, USNO, CAHA 2.2-m, ESO VLT+FORS1, 
ESO 2.2-m+WFI in the optical bands, and the UKIRT obser- 
vations in the infrared bands. The magnitudes were obtained 
from PSF photometry of the OT using DAOPHOT II, except 
for the VLT observation of 2000 March 6.39 UT, which has 
been derived from the March 6 VLT spectra (see Table ph. 



Telescope 


date 


seeing 


filter 


exp. time 


mag 




(2000. UT) 


ll'C ^cc 




(sec) 




NOT 


IVldl J.Jt 


2 3 




900 


on OQ+O 04 


NOT 


Mar 3.17 


2.2 


R 


900 


20.15±0.04 


NOT 


Mar 3. 19 


2.1 


R 


900 


20. 11 ±0.04 


NOT 


Mar 3.20 


2.2 


B 


900 


2 1.05 ±0.04 


NOT 


Mar 3.21 


2.3 


R 


900 


20.14±0.04 


NOT 


Mar 3.22 


2.3 


B 


900 


21.02±0.04 


NOT 


Mar 3.24 


2.3 


U 


900 


20.52±0.08 


NOT 


Mar 3.25 


2.2 


R 


900 


20.16±0.04 


NOT 


Mar 3.26 


2.0 


I 


900 


19.60±0.07 


NOT 


Mar 3.27 


2.3 


I 


300 


19.59±0.06 


CAHA 2.2m 


Mar 4. 14 


1.4 


I 


1440 


20.10±0.10 


USNO 


Mar 4.39 


2.0 


I 


480 


20.04±0.11 


USNO 


Mar 4.40 


2.1 


V 


720 


21.05±0.07 


USNO 


Mar 4.41 


2.0 


B 


1200 


21.35±0.05 


USNO 


Mar 4.42 


1.8 


R 


600 


20.61 ±0.06 


USNO 


Mar 4.48 


1.6 


R 


1200 


20.58±0.03 


USNO 


Mar 4.49 


1.6 


R 


1200 


20.54±0.04 


USNO 


Mar 4.50 


1.6 


R 


1200 


20.60±0.04 


USNO 


Mar 4.52 


1.6 


B 


1200 


21.42±0.04 


USNO 


Mar 4.53 


1.6 


V 


1200 


20.98±0.04 


USNO 


Mar 4.43 


2.0 


U 


1200 


20.82±0.11 


UKIRT 


Mar 4.51 


0.7 


J 


810 


19.19±0.04 


UKIRT 


Mar 4.53 


0.6 


H 


810 


18.44±0.04 


UKIRT 


Mar 4.55 


0.5 


K 


810 


17.62±0.05 


CAHA 2.2m 


Mar 5.23 


1.5 


I 


960 


20.24±0.12 


ESO VLT 


Mar 5.39 


0.9 


R 


10 


20.61±0.05 


ESO VLT 


Mar 6.39 


1.1 


R 




21.43±0.26 


NOT 


Mar 7.22 


1.4 


R 


300 


21.59±0.07 


NOT 


Mar 7.24 


1.4 


U 


3600 


21.92±0.07 


CAHA 2.2m 


Mar 7.23 


1.1 


I 


960 


21.13±0.15 


NOT 


Mar 8.18 


1.4 


R 


1200 


21.80±0.05 


NOT 


Mar 8.20 


1.5 


I 


900 


21.39±0.09 


NOT 


Mar 8.21 


1.6 


B 


1100 


22.89±0.08 


NOT 


Mar 8.25 


1.9 


U 


3700 


22.49±0.12 


NOT 


Mar 9.15 


1.4 


R 


2000 


22.11±0.15 


NOT 


Mar 9.20 


1.6 


U 


5400 


22.69±0.20 


NOT 


Mar 9.24 


1.4 


B 


1200 


23.01±0.11 


NOT 


Mar 9.26 


1.2 


I 


900 


21.72±0.15 


CAHA 2.2m 


Mar 10.05 


1.3 


I 


960 


21.85±0.2() 


USNO 


Mar 10.40 


3.0 


B 


1800 


23.40±0.3() 


CAHA 2.2m 


Mar 11.21 


1.2 


I 


1440 


22.63±0.27 


ESO 2.2m 


Mar 11.39 


2.4 


R 


3200 


23.12±0.18 


USNO 


Mar 12.44 


1.9 


R 


4800 


23.10±0.22 



4.1. Optical data 

To avoid contamination from the nearby star A (located at a 
separation of 6" west and 1" south of the OT in Fig. ^J), we 
measured the magnitude of the OT relative to stars in the field 
by performing Point Spread Function (PSF) photometry, using 



DAOPHOT II (Stetson 1987, 1997). There are several bright 



4. Photometry 



and unsaturated stars in the field from which a good PSF could 
be determined. 

For the data presented here, there is no indication of a con- 
tribution from a hos t ga laxy to the emission at the position 
of the OT (see Sect. 1J1 for a discussion of the host galaxy). 
Hence, extended emission from a faint galaxy at the position of 
the OT will not affect the PSF photometry appreciably (much 
less than observational errors). The quality of the PSF photo- 
metry was checked by subtracting the PSFs from the images 
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of star A and the OT. In all frames the residuals are consistent 
with being shot-noise. 

To avoid errors due to colour terms or colour differences in 
our photometry (the conditions during most of the observations 
at the NOT were possibly non-photometric due to increasing 
amounts of Saharan dust in the atmosphere over the telescope at 
the time of observations), the magnitudes of the OT for all our 
optical photometry were calibrated relative to stars of similar 
colours in the field. 

The photometric standard UBVRcIc calibration of the 
field was performed at the USNOFS 1.0-m telescope and is 
available in Henden (200C). This calibration has an estimated 
zero-point uncertainty of 2 percent, which is well below the 
errors in the relative magnitudes. The results of the PSF pho- 
tometry are presented in Table |l[ The 2000 March 6.39 VLT R- 
data point has been derived from the March 6 combined VLT- 
spectrum (Table ||). 

Based on this photometric calibration we conclude that 
star A showed no sign of variability within observational er- 
rors throughout our observations, and that it had the fol- 
lowing magnitudes: U = 20.427±0.133, B = 19.837±0.030, 
V = 18.767±0.018, R = 18.084±0.043, and I = 17.526±0.044. 

4.2. Near-infrared data 

The UKIRT images were processed using the ORAC imag- 
ing data reduction routines developed for UKIRT (Bridger et 
al. 2000 ). The J, H and K magnitudes of the OT were then mea- 
sured from the UKIRT data as follows. First we measured the 
magnitude of the OT relative to star A using DAOPHOT II PSF 
photometry as described above. Then, in order to transform 
this magnitude to the standard UKIRT system, we performed 
aperture photometry in an aperture with a diameter of 2'.'7 on 
calibration images obtained of the standard stars S868-G and 
p389-d from the list of UKIRT faint standards^] and on star A. 
The estimated error in the zero-point is about 0.05 in each of 
J, H and K. We have assumed negligible extinction difference 
between standard and program field. The results of these mea- 
surements are presented in Table [j] 

5. Spectroscopy 

Spectroscopic observations were carried out on 2000 March 5 
and 6 UT with VLT-Antu equipped with FORS1 (for details 
see the observing log in Table g). We used the GRIS_300V+10 
grism and the GG375 order separation filter, which provide 
a spectral coverage from 3600 A to 8220 A and a disper- 
sion of 2.64 A/pixel. The effective exposure time was 800 s 
on March 5.39 UT and 1200 s on March 6.38 UT. Standard 
procedures were used for bias and flat field correction, and 
the opti mal ex traction procedure for faint spectra (described in 
M0ller ( 2000 )) was used to extract one dimensional spectra. 

The position angle of the long slit was chosen such that 
both the OT and star A were centred onto the slit. From the 
magnitude of star A we calibrated the flux of the trace of the 



Table 2. Spectroscopy of the afterglow of GRB 00030 1C, ob- 
tained with VLT+FORS1 on 2000 March 5 and 6 UT. 



Date 


Exposure time 


Airmass 


Slit width 


(2000 UT) 


(s) 




(") 


Mar 5.39 


400 


1.732-1.745 


1.3 


Mar 5.40 


400 


1.719-1.729 


1.3 


Mar 6.38 


600 


1.772-1.801 


1.0 


Mar 6.39 


600 


1.745-1.767 


1.0 






Fig. 4. Combined VLT+FORS1 spectrum of GRB 00030 1C 
from 2000 March 5+6 UT. The spectrum is normalized to 1 
in the continuum. The atmospheric absorption bands and resid- 
uals from strong sky-lines are marked, as well as the 4 ab- 
sorption lines listed in Table [}| The spectrum is binned to 7 A 
pixels, and the lower curve shows the noise (per pixel). 

optical counterpart. The spectral flux calibration derived for the 
OT on March 5 is consistent with the optical photometry dis- 
played in Table [j]. We derive a value for the spectral index of 
13 = -1.15 ± 0.26 on March 5 and [3 = -1.43 ± 0.28 on 
March 6, corrected for interstellar extinction of E(B — V) = 
0.053 ± 0.020, using the dust maps of Schlegel et al. ( |1998| ). 



|http: // www . jach . hawaii . edu/ JACpublic/UKIRT/| 



astronomy /calib/f aint_stds . html 



6. Results 

6. 1. GRB absorption lines and redshift 
determination 

The combined spectrum, reproduced in Fig 3, has a resolution 
of 14 A FWHM, and signal-to-noise (S/N) in the range 15-30 
per resolution element redwards of 4000 A. From 4000 A to 
3600 A the S/N drops rapidly. 

Due to the poor resolution, only very strong absorption 
lines can be detected individually. In Table || we list the only 
four absorption features which were detected at a S/N in excess 
of 4.5. Two of the features were found bluewards of 4000 A, 
and were initially ignored. The line at 4712 A is broader than 
the resolution profile, and we tentatively identified it as a pos- 
sible C IV absorption complex with redshifts in the range 2.038 
to 2.042. With this identification, the other three features would 
fit the proposed identifications given in Table ||. Note, however, 
that the Si II 1260 line is far too strong and wide to be a sin- 
gle line, and we hence assume that it is a blended feature. The 
Fe II 2600 line is strong and narrow, but also this line seems 
excessively strong given the lack of other strong Fe II lines. 
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Table 3. Absorption features detected at S/N>4.5 in the spec- 
tram of GRB 00030 1C on 2000 March 5+6 UT. W obs indicates 
the observer frame equivalent width of the features. 



Line 


Identification 


Ayac 


w obs 


<t(W) 




(A) 


(A) 






1 


Lya 


3693.56 


67.07 


12.09 


2 


Sill 1260 


3843.09 


18.99 


2.03 


3 


Civ 1549 


4712.23 


3.76 


0.82 


4 


Fe II 2600 


7909.03 


6.29 


1.21 



In order to provide a more strict test of our proposed iden- 
tification, and to obtain an accurate value for the redshift, we 
proceeded as follows. First we shifted and stacked pieces of 
the spectra where we would expect common low ionization 
lines. We selected the singly ionized species of Si, C and Fe, 
all of which are known as strong absorbers in quasar absorp- 
tion line systems. In total our spectrum covers positions of the 
lines Fe II 1608, Fe II 2344, Fe II 2374, Fe II 2382, Fe II 2586, 
Fell 2600, Sill 1260, Sill 1304, Sill 1526, and C II 1334. 
Si II 1260 (at 3832 A) was in the very low S/N part of spectrum, 
and almost certainly blended, so it was not included. Treating 
each ion separately, a weighted mean absorption feature was 
calculated using the oscillator strength of each line as statisti- 
cal weight. 

The regions of the co-added spectra for each of the ions 
Fe II, Si II and C II, transformed into redshift space, is shown in 
the left panel of Fig. |[ A combined "Low Ionization" absorp- 
tion feature (bottom of left panel of Fig. |[) was obtained by co- 
addition of the three sets of features, but using the number of 
lines as statistical weights. The redshift range searched for low 
ionization absorption systems by this method was 1.95-2.14 
and no other candidate systems were found in this range. For 
comparison of the redshifts, we plot in the right panel of Fig. |] 
again the combined "Low Ionization" absorption feature to- 
gether with the C IV absorption trough. The bottom panel here 
shows the combination of all the lines. Given the significance 
of this combined "Metal Absorption Feature", we conclude that 
the tentative identification of this system is confirmed. 

It is commonly seen in quasar absorption line systems that 
the low ionization species, tracing the cold dense gas, have 
a more well-defined redshift than the high ionization species. 
Hence, we shall adopt the redshift z abs = 2.0404 ± 0.0008, 
measured from the combined Low Ionization feature, as the 
systemic redshift. It is seen from the bottom right panel of 
Fig. H that inclusion of C IV would result in a slightly higher 
redshift. Note that z^bs = 2.0404 ± 0.0008 is consistent with 
the redshift z=1.95±0.1 based on the Lyman break (Smette et 
al. 2000; Feng et al. 2000), but is significantly hi gher th an the 
value z=2.0335±0.0003 reported by Castro et al. (§00(J. 

The oscillator strength weighted mean observed equiva- 
lent width of the Fe II lines is 2.56 A, which is strong enough 
that by comparison to known quasar absorbers one would ex- 
pect this to likely have a column density of neutral Hydrogen 
in excess of 2 x 10 20 cm~ 2 . Such absorbers are known as 
Damped Lya Absorbers (DLAs), and hold a special interest 
because of the large amounts of cold gas locked up in those 
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Fig. 5. Left: Upper panel: Oscillator strength weighted mean 
absorption lines of Fe II, Si II and C II, lined up in redshift 
space for easy comparison (see text for details). Lower panel: 
The combined "Low Ionization" absorption feature (weighted 
mean of the three sections of spectrum shown above). Right: 
Upper panel: The combined "Low Ionization" absorption fea- 
ture lined up with the C IV trough. Lower panel: Combined 
"Metal Absorption Feature". All the sections of spectra have 
been normalized to 1 in the continuum, but shifted along the 
abscissa by suitable off-sets. Note that none of the figures start 
at 0, instead the scale on the abscissa provides the proper refer- 
ence. 



objects (Wolfe et al. |1995| ; Storrie-Lombardi et al. |1997| ). It 
is commonly assumed that the DLAs are the progenitors of 
present day disk galaxies, but they have proven extremely dif- 
ficult to identify (see e.g. M0ller & Warren 1993; Kulkarni 
et al. 2000b). Observational evidence has been accumulating 
(M0ller & Warren 1998; Fynbo et al. |1999| ) which suggests 
that a likely reason why DLA galaxies are so hard to identify 
is their small gas cross-sections and faint magnitudes, caus- 
ing them to stay hidden under the point spread functions of 
the bright quasars. A GRB selected DLA galaxy sample would 
not be hampered by this problem once the OT has faded suf- 
ficiently, and could as such help greatly in understanding the 
nature of DLA galaxies. We shall therefore now briefly con- 
sider the low S/N part of the GRB spectrum below 4000 A, to 
investigate if any information concerning the H I column den- 
sity can be extracted. 
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Redshift 

Fig. 6. Upper panel: The "Combined Metals" absorption fea- 
ture from Fig. |] is shown here for comparison. Lower panel: 
Spectral region around Lya (for comparison to the metal ab- 
sorption also plotted in redshift space). Note the very sharp 
onset of absorption well above the expected redshift. This is 
consistent with a very broad Lya absorption line as detailed in 
the text. The lower curve shows the noise per bin (bins of 0.003 
in redshift). 



In Fig. H (lower panel) we have plotted the spectral region 
around Lya, and for comparison of redshifts, the "Combined 
Metals" feature (upper panel). Also plotted on the lower panel 
is the noise per bin (for redshift binsize = 0.003). It is clearly 
seen that the spectrum drops steeply before the expected cen- 
tral position of the Lya line, and well before the S/N drops 
below detection. One likely explanation for this is the pre- 
sence of a very broad Lya absorption line. To quantify this 
we have modelled several Lya absorption lines, all at redshift 
2.0404, and calculated the \ 2 of their fit to the data in the range 
3700 A to 3750 A. For N(H i) = the X 2 per degree-of- 
freedom DOF is 6.46 which confirms that an absorption fea- 
ture is indeed present. The formal % 2 minimum is found at 
N(H i) = 1.5 X 10 21 cm- 2 (x 2 per DOF = 0.86), but any value 
within a factor 3 of this is acceptable. 

It should be recalled that the above estimate, in a strict 
sense, only applies in the case the OT lies "behind" the absorb- 
ing cloud. In case the OT is in fact embedded "inside" a DLA 
cloud, resonant scattering of Lya photons may alter the profile 
of the absorption line somewhat. In the present case this is a 
detail which the quality of our data will not allow us to discern. 



6.2. The multi-wavelength spectrum around March 
4.5 UT 

The wide wavelength coverage obtained around 2000 March 
4.5 UT allows us to construct the multi-wavelength spectrum of 
the afterglow, by using our USNO and UKIRT data. Mos t fire 
ball models (e.g. Sari et al 



1998 


! Piran 1999 


l Meszaros 


1999 



and references therein) and observations of previous afterglows 
suggest a power-law Spectral Energy Distribution (SED). 
However, a glob al fit to the broadband SED of GRB 00030 1C 
by Berger et al. ( 2000 ) demonstrates that the mm-optical range 
cannot be described by a single power-law. Thus, we have only 
considered wavelengths shorter than IR when fitting the SEDs 
presented in this section. This gives eight measurements in the 
period March 4.39^.55 UT, namely UBVRIJHK. The eight 
measurements, plotted against \ogv together with the VLT 
spectrum in Fig. [7| allow us to constrain the SED of the OT 
at this epoch. 

As extinction is highly wavelength dependent, a progres- 
sive deviation from a pure power-law fit can be explained as 
being due to dust extinction in the GRB host galaxy. In or- 
der to test this possibility, we first fit a power-law to the data 
in the NIR-optical range with the addition of an intrinsic ex- 
tinction law, i.e., using the expression F v oc v& x 10(~°- 4j4 "), 
where A v is the rest frame extinction in magnitudes at the rest 
frame frequency v. Four extinction laws have been applied in 
order to establish a relatio nship between A„ and v. Following 
Rhoads & Fruchter (2001 1, we have fitted the extinction curves 



for t he Ga laxy and Magellanic Clouds published by Pei ( |1992[ ). 
Pei ( 1992 ) provides extinction laws for the Milky Way (MW), 
the Large Magellanic Cloud (LMC), and the Small Magellanic 
Cloud (SMC) based on their different proportions in the dust- 
to-gas ratio (1:1/5:1/10) and in the abundance of heavy ele- 
ments (1:1/3:1/8). The most significant difference is the se- 
quential change in the strength of the 2175 A extinction feature, 
prominent for the MW, moderate in the LMC, and nonexistent 
for the SMC extinction curve. It is important to note that for 
the redshift of GRB 00030 1C this extinction feature falls in the 
observed R-band. So, the presence of this feature would result 
in a clear decrease in the R-band flux compared to the I and 
V-bands. The MW extinction curve requires an equal amount 
of graphite and silicate grains, while the SMC extinction curve 
can be explained by silicate grains only, with the LMC extinc- 
tion law as an intermediate stage. 

We have applied these four extinction laws to our data in 
order to infer qualitative information about the dust-to-gas ra- 
tio, the abundance of heavy elements and the composition of 
the dust in the host galaxy of GRB 00030 1C. We leave A„ as a 
free parameter, so fitting a function like F u oc v& x 10^ OAA ^ 
allows us to determine (3 and A v simultaneously. The values of 
X 2 are displayed in Table |] for each case. 

A pure power-law fit (F v oc ir~) to the eight data points 
(after correction for Galactic extinction) leads to a reduced % 2 
of 1.69 (see Table Q), making an acceptable description of the 
NIR-optical range of the SED. However, the fit can be im- 
proved if a modest amount of extinction is introduced. This 
is because the SED is slightly bending down towards higher 
frequencies (see Fig. 0). The eight data points show that there 
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Fig. 7. Spectral energy distribution of the OTon 2000 March 4.39-4.55 UT (log v— log F„). The frequency scale has been redshift 
corrected to rest frame values assuming z = 2.0404 (see Sect. 6. 1 ). The data points include our optical (UBVRI) and near infrared 
(JHK) data (see Table [l]) from March 4.5 UT. In the plots, a range of extinction models have been fitted to the data, cf. the text 
and Table || Upper panel: UBVRIJHK photometry with the VLT spectrum (dotted line) of March 5.4 included (the offset from 
the photometry is due to the fading of the OT between March 4.5 and 5.4). The solid line represent the fitted SED when the 
emperical SMC extinction law is applied. A fit (longdashed) assuming no extinction of the host is shown for comparison. Lower 
panel: here is showed the effects on the LMC and MW-like SED fits when reasonable values of (3 and Ay are considered. The 
solid line, as in the upper panel, shows the SED when a SMC-like extinction law is fitted. Then, the values obtained of j3 and Ay 
{(3 = 0.7, Ay = 0.09) are applied to plot power-law afterglow SEDs with MW (long dashed) and LMC-like (dashed) extinction. 
Circles: USNO UBVRI-photometry, Triangles: UKIRT JHK-photometry. 



is not any presence of a redshifted 2175 A absorption bump in 
the R-band at aU. In short, the near-IR SED of GRB 00030 1C 
can be described as a curved power-law but with no broad ab- 
sorption features. 



As expected by the lack of the absorption bump in the R- 
band, the MW and the LMC extinction laws are completely in- 
consistent with our data. In fact, both fits imply an unphysical 
negative extinction (see Table 0). This is because the R-band 
flux is slightly over the linear interpolation between the I-band 
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Table 4. Result of fitting different extinction law models to ob- 
servations of the afterglow on 2000 March 4.5 UT. The ex- 
tinction curves include the models by Pei (1992) for the Milky 
Way (MW), the Large Magellanic Cloud (LMC) and the Small 
Magellanic Cloud (SMC) (see Fig. and text for further de- 
tails). 





X^/DOF 


13 


A v 
mag 


No extinction 


1.69 


-0.90± 0.02 





MW, Pei (1992) 






< 


LMC, Pei (1992) 






< 


SMC, Pei (1992) 


0.91 


-0.70± 0.09 


0.09 ± 0.04 



Table 5. Best fits for the colours of the OT from 2000 March 3 
to 1 1 UT, assuming an achromatic evolution. V— R is for March 
4.4 only. P(x 2 ) is the probability to obtain a lower value of 
X 2 /DOF for the given model (constant colour). Colours are not 
corrected for galactic or intrinsic extinction. See also Fig. H. 



Colour 


Value 


a 


X 2 /DOF 


P(x 2 ) 


U-B 


-0.49 


0.06 


0.41 


0.26 


B-R 


0.91 


0.03 


1.67 


0.83 


R-I 


0.49 


0.04 


0.51 


0.20 


B I 


1.43 


0.05 


0.33 


0.14 


U-I 


0.92 


0.07 


0.64 


0.36 


U-R 


0.40 


0.05 


1.35 


0.75 


V-R 


0.44 


0.05 







and V-band fluxes (in a Log-Log space), and both of these two 
extinction laws fit the 2175 A bump as an emission feature in- 
stead of an absorption bump. To illustrate the problem with the 
MW and LMC extinction curves, we have, in the lower panel 
of Fig. plotted the effect of having the MW and LMC extinc- 
tion laws with the parameters derived for the SMC extinction 
law (/3=-0.7, Ay=0.09). As seen, the shapes of both SEDs are 
incompatible with our UBVRIJHK measurements. 

The quality of the MW, LMC and SMC and unextincted 
SEDs can also be compared checking the flux pr edicte d at 
250 x (1+z) GHz (rest-frame), where Berger et al. ( |2000| ) re- 
ports a flux of 2.1 ± 0.3 mJy at 250 GHz on March 4.29 UT. 
Extrapolating the four extinction curves, we obtain the follow- 
ing fluxes in increasing order; 5.0 mJy (SMC), 18.1 mJy (No 
extinction), 33.9 mJy (MW) and 51.6 mJy (LMC). As with the 
NIR-optical range, the SMC extinction provides the most rea- 
sonable results. The actual measured flux in mm is below the 
value predicted by the mm-optical extrapolation, because the 
pure power-law assumption is not correct in the mm-NIR spec- 
tral range and an additional curvature effec t is present in the 
SED, as demonstrated by Berger et al. (2000 1 (see their Fig. 1). 
Thus, the value of Ay = 0.09 ± 0.04 for the SMC-fit given in 
Table ^| should be taken as a good indication of the real extinc- 
tion, although strictly speaking it is just an upper limit. 

In conclusion, the featureless SMC extinction law provides 
the best fit to our data, improving the quality of the fit ob- 
tained for an unextincted afterglow (see Table ||). It is inter- 
esting to note the dramatic dependence of the quality of fit on 
the existence of the 2175 A absorption bump. Extinction laws 
with high moderate dust-to-gas ratios that produce such an ab- 
sorption feature do not provide good fits to our data points. 
Therefore, the spectral energy distribution of GRB 00030 1C 
supports a scenario where the host is in the early stages of 
chemical enrichment. 

The power-law + extinction fits to the SED in the NIR- 
optical range allow us to predict the UV flux at the spec- 
tral range of MAMA-HST when the UV spectrum was ob- 
tained at 2000 March 6.375 UT (Smette et al. |2000b assum- 



ing that the shape of the SED has not changed between the 
two epochs. This assumption is supported by the imaging data 
(Sect. 6.3). First, we consider the best fit to the SED at 2000 
March 4.39-4.55 UT and calculate the flux at March 4.39-4.55 
UT at 3000 A. Then, making use of the light curve models 



(presented in Table |5|), we estimate the value of the flux at 
3000 A for March 6.375 UT. The predicted flux ranges from 5.9 

-2, 



x 10 18 erg cm * s 



- 1 A" 1 to 7.9 xl0~ 18 erg cm- 2 s" 1 A" 1 , 



depending on the light curve model. A final analysis (Smette 
et al. |2001 1) of the MAMA-HST data revealed a flux of ~ 
erg cm -1 consistent with our extrapolation. 



et al. pOOl 

7.3±?;| 10 



6.3. Evolution of the spectral energy distribution 

From our UBRI photometric data (presented in Table |]) we 
have multi-band optical coverage from 2 to 10 days after the 
burst-trigger (on March 1.41 UT). When analysing the colours, 
we find that the simplest reliable fit is for constant colours. 
Thus we find no evidence for optical chromatic evolution for 
the afterglow during the period of observations (see Fig. |]). 
For these constant fits, we obtain the values presented in 
Table ||. These values are not corrected for Galactic or intrinsic 
extinction. 



6.4. The light curve 

According to the simple fireball model the optical afterglow 
should follow a power-law decay, F v oc v^t a (Sarietal. 1998). 
However, a single power-law is excluded at more than the 
99.9% confidence level. The parameters for this power-law 
are given in Table ||. The photometry suggests that the opti- 
cal afterglow follows a shallow power-law decay for the first 
few days and then steepens. This behaviour has been seen pre- 



viously in GRB 9805 19 (Jaunsen et al. |2001j), GRB 99012 3 
(Kulkarni et al. p999b, GRB 990510 (Harrison et al. p999b, 



GRB 991208 (Castro-Tirado et al. 2001) and GRB 000926 



(Fynbo et al. 2001b and is predicted by many models for 
gamma-ray bursts (see below). Sagar et al. (2000) report that 
there are seven components to the R-band light curve. Here 
we are primarily interested in the overall structure of the light 
curve, not the structure at small time scales. Therefore, we fit a 
broken power-law of the form, 



fu(t) = 



f»(tb)U) a \ if t <t b 



(i) 
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Fig. 8. Fitting the colours from the optical photometry from 
2000 March 3 to March 1 1 UT, assuming achromatic evolu- 
tion. The solid line indicates the best fit, with the dashed lines 
representing the error-range. The photometry is from Table [j]. 
The results of the fitting is presented in Table |[ 

to the UBRI data presented in Table [j]. The flux, in /iJy, at time 
t days after the burst is denoted by f v (t). The time of the break 
in the decay is denoted fa- The slope before the break is ot\, 
and the slope after the break is ot-x. The flux at the time of the 
break is f v {fa). We used CERN's MINUIT function minimiza- 
tion package, and a chi-square minimization scheme, to simul- 
taneously solve for the four free parameters (ai, a 2 , fa, and 
f v (fa)) and their formal l-a errors in the fit for each parameter. 

The data was corrected for Galactic reddening and extinc- 
tion before the fits were made. No corrections were made for 
reddening or extinction in the host galaxy. The photometry was 
transformed to the R band using the colours given in Table | 
and then converted to units of flux using a photometric zero 
point of f v % o = 3.02 x 10~ 20 erg cm -1 s -1 Hz T 1 (Fukugita 
et al. |1995| ). 

The best fit to the combined UBRI photometry is listed in 
Table |(| and shown in Fig. |[ To test the sensitivity of the results 
to the fitting function, we a lso fit our data with the smooth func- 
tion used by Stanek et al. ( |l999t their Eq. 1) on GRB 990510. 
The results are given in Table |6j The broken power-law gives 
the smallest chi-square value, and the errors in the individ- 
ual parameters are smaller for the broken power-law fit than 
they are for the smooth function. The correlation coefficient 



Table 6. The parameters of the best-fitting functions to the 
optical decay of GRB 00030 1C. The number of degrees-of- 
freedom (DOF) in each fit is the number of data points minus 
the number of parameters. The number of parameters is the 
number of free parameters in each model plus the number of 
colours that the data was adjusted with, in order to bring it into 
the R band. 





Single Power Law 


Broken Power Law 


Smooth Function 


a 


-1.22 ± 0.07 






ai 




-0.72 ±0.06 


-0.55 ± 0.14 


C«2 




-2.29 ± 0.17 


-3.53 ± 0.47 


h 





4.39 ± 0.26 


5.99 ± 0.79 


/,(*„) 


70.46 ± 0.56 


16.23 ± 1.13 


8.24 ± 2.05 


Ro(t b ) 


19.08 ± 0.01 


20.67 ± 0.08 


21.41 ± 0.27 


X 2 /DOF 


4.499 


1.023 


1.090 


DOF 


33 


31 


31 



between fa and a± is —0.39 and the coefficient between fa and 
«2 is —0.84. The broken power-law fit is consistent with the 
data at the 43% confidence level. Even though we, from the 
theory of fireballs, would expect that the light-curve evolution 
is a smooth function, we find in the case of GRB 00030 1C, in 
agreement with Berger et al. ( |2000| ), that the broken power-law 
provides the most reliable fit. Additionally, a broken power-law 
provides the most reliable metod of determining the time when 
the decay of the light has steepened, and thus is a useful way 
of parameterising the data. 

We combined all our UBRI-data (Table [p) with those from 
the literature. Fig. shows the best-fitting broken powe r-law 



for all of the UBRI data in the literature (Sagar et al. |2000 
and references therein) and Table |l[ This data was shifted to 
the R band in the manner described above. The parameters of 
the fit are shown in Fig. [l(] and are not significantly different 
from the parameters of the broken power-law that was fit to 
our data (see Table ^). The conspicuous short-term be haviou r 
of the light -curve has been detaile d by M asetti et al. ( 2000 ), 
Sagar et al. (|200C| ) and Berger et al. ( |2000[ ). Garnavich, Loeb & 
Stanek ( 2000 ) find that the variation of the lightcurve can be in- 
terpreted as a microlensing event, peaking about 3.5 days after 
the burst, superposed on a power-law broken at fa = 7.6 days. 
This superposed event peaks at a more sparsely sampled period 
in our data, coinciding partly with where we identify the break. 
Thus it is not possible, from our data, to further constrain the 
existence of such an event. We choose here to work only with 
the data presented in Table |l] as it represents a consistently de- 
rived set. 



7. Discussion 

7. 1. Interpretation of the light curve 

The fit to the multi-colour light curves shows a break at fa = 
4.39 ± 0.26 days, with the light curve steepening from u\ = 
-0.72 ±0.06 to a 2 = -2.29 ± 0.17, i.e., by Aa = ai - a 2 = 
1.57 ± 0.18. A broken light curve can arise in a number of 
circumstances: i) If the frequency separating fast cooling elec- 
trons from slow cooling ones moves through the optical at fa, 
the resulting light curve would steepen by Aa ~ 0.25 (Sari 
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Fig. 9. The upper panel shows the UBRI photometry from 
Table [j] and the best-fitting broken power-law fit to the data. 
The photometry was offset to the R band using the colours 
given in Sect. 5.3 (see the text). The horizontal line at the lo- 
cation of the break shows the l-a uncertainty in the time of 
the break. The lower panel shows the residuals of the fit in the 
sense R b s — Rfit- The uncertainties in the residuals are the 
uncertainties in the photometry. 
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Fig. 10. This Figure is analogous to Fig. but instead 
shows the data from Table [l] (filled circles) together with all 
the UBR I pho tometry presented in the literature (Rhoads & 
Fruchter 2001 and references therein) (open circles) and with 
the best-fitting broken power-law fit to this data. 



et al. 1998| ). ii) The light curve may also steepen if a spherical 
fireball slows down to a non-relativistic expansion (Dai & Lu 



1999), resulting in Aa = — (ai+3/5) = 0.12 for our value of 



ai. Hi) If the outflow is collimated with a fixed opening angle, 
the break in the light curve occurs when the relativistic beaming 
of the synchrotron radiation becomes wider than the jet open- 
ing angle (Meszaros & Rees 1999). In this case the break is 
a geometrical effect and the steepening is Aa = 3/4. iv) If 
the afterglow arises in a sideways expanding jet, the steepening 



will be Aa = (1 - ax/3) = 1.24 (Rhoads |T999| ) for our value 
of a±. The above estimates all assume a constant mean density 
distribution of the ambient medium. We note that collimated 
outflows in general result in faster decaying light curves than 
the spherically symmetric ones. If the mean density distribution 
is not constant, e.g. it has a stellar wind density profile, the light 
curves also deca y faste r, but the break will be less pronounced 
(Panaitescu et al. |1998| ). 

Based on the light-curve properties alone, the model that 
best fits the observations is that of a sideways expanding jet in 
an ambient medium with a constant mean density distribution. 
In that interpretation, the observed light curve indices imply 
an electron energy distribution index of p = 2.13 ± 0.09 that 
results in a theoretical spectral index of (3 = — (p — l)/2 = 
— 0.56±0.05. This is in agreement with the spectral index of 
(3 = — 0.70±0. 09 inferred from our spectroscopic observations 
when correcting for extinction in the host galaxy (Sect. 6.2), 
independently strengthening the described model for the after- 
glow. 

With a combined fluence in the 25 - 1000 KeV range of 
4 x 10~ 6 erg cm~ 2 , and a redshift of z = 2.0404, the isotropic 
energy release of GRB 0030 1C is E = 4.6 x 10 52 erg. 
Following Rhoads dl999b , the energy estimate and the light 
curve break time, tb — 4.39 ± 0.26 days, implies a jet opening 
angle, at that time, of 9 ~ IS ??. 1 / 8 , where n is the number den- 
sity of the ambient medium (in units of cm -3 ), and the break 
is assumed to occur when the opening angle equals the bulk 
Lorentz factor, 6 = T. 

This interpretation is similar to that of GRB 990510 which 
was almost 5 times more energetic, but had a jet opening angle 
of approximately 5°, leading to an earlier break in the light 
curve. The best fit to GRB 990510 was a smooth function (e.g. 
Stanek et al. |1999|; Harrison et al. |1999|; Holland et al. |200C(), as 



compared to a broken power-law in the case of GRB 00030 1C. 



7.2. The host galaxy 

The host galaxy appears to be very faint. There is no evidence 
for any extended emission from a host galaxy in any of the data 
presented in this paper. Deep images obtained with HST+STIS 
about 1 month after the GRB indicate that any host galaxy 
must be R>27.8±0.25 (Fruchter et al. |2000| ). Hence, we can 
safely conclude that the host galaxy of GRB 00030 1C is very 
faint compared to other known populations of galaxies at high 
redshifts. From fitting extinction curves to our photometry, we 
have found evidence for some extinction (Ay ~ 0.1) in the 
host galaxy. The Ay derived from the best fit corresponds to an 
absorption at rest-frame 1500 A of about 0.4 magnitudes. 
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For comparison, Lyman Break Galaxies (LBGs) at slightly 
higher redshifts (z ~ 2.5-3.5) on average have values of ex- 
tinction at rest frame wavelength 1500 A of approximately 
1.7 magnitudes, and, in rare cases up to 5 magnitudes (Steidel 
et al. 1999| ). Hence, the faint optical appearance of the host 
galaxy relative to the star-forming LBGs at high redshift is 
most likely not imposed by massive extinction, but is rather 
due to a lower overall star formation rate of the host galaxy. 

As the host galaxy furthermore has a very high H I column 
density, log(NH l)=21.2±0.5 as derived from the Lya absorp- 
tion f eature and supported by the strong Lyman break (Smette 
et al. 2000| ), it is interesting to compare with the population 
of galax ies identified as Damped Lya Absorbers (Wolfe et 
al. 1986| ) in the spectra of background QSOs. These galaxies 
have Hi column densities higher than log(NHl)=20.3. Based 
on the luminosity function of LB Gs an d the typical impact pa- 
rameters of DLAs, Fynbo et al. ( 1999| ) show that the majority 
of DLAs at z = 3 must be fainter than the current flux limit for 
LBGs of R=25.5 and that there hence is a very abundant popu- 
lation of galaxies fainter than the LBG flux li mit. A s imilar con- 
clusion has been reached by Haehnelt et al. ( 2000 i. The dust- 
to-gas ratio towards the line-of-sight of GRB 00030 1C gives 
a value of A y /N(Hi) < 0.1/1.6 x 10 21 cm 2 = 0.6 x 10~ 22 
cm -2 . This upper limit is, within errors, consistent with the ex- 
pected Ay/N(H i) for DLAs. The corresponding value for the 



Milky Way is 2 x 10" 22 cm" 2 (Allen 2000) 



It is still uncertain what fraction of the integrated star for- 
mation rate at high redshift is accounted for by the LBGs and 
what fraction has to be accounted for by galaxies further down 
the luminosity function. The relative occurrence of GRBs in a 
given population of galaxies is expected to be proportional to 
its rel ative c ontribution to the to tal star format ion ra te (Totani 
et al. |l997j Wijers et al. |l998t Mao & Mo |l999fc Blain & 



Natarajan pOOOP . However, so far only one GRB (GRB 971214 
at z = 3.418) is confirmed to have occurred in a galaxy simi- 
lar to the faint members of the LBGs selected in ground based 
surveys (Odewahn et al. |1998| ). The fact that GRB 000301C oc- 
curred in an intrinsically very faint galaxy and that most GRBs 
with identified OTs have occurred in L* or sub-L* galaxies, 
suggest that a large fraction of total star formation at high red- 
shift occurs in a population of galaxies that is further down 
the luminosity function than the bright LBGs found in ground 
based surveys and that is likely to have a large overlap with the 
DLAs. 

8. Conclusion 

GRB 00030 1C is so far the GRB of shortest duration, for which 
a counterpart has been detected. The high-energy properties of 
the burst are consistent with membership of the short-duration 
class of GRBs, though GRB 00030 1C could belong to the pro- 
posed intermediate class of GRBs or the extreme short end of 
the distribution of long-duration GRBs. Our VLT-spectra show 
that GRB 000301C occurred at a redshift of 2.0404±0.0008. 
The light curve of the optical transient is well-fitted by a 
broken power-law and it is consistent with being achromatic. 
From the light-curve properties we find that the best model for 
GRB 00030 1C is that of a sideways expanding jet in an ambi- 



ent medium of constant density. This interpretation is further 
supported by the achromatic light-curve evolution, and by the 
agreement between the theoretically predicted and observation- 
ally derived spectral indices. The spectral energy distribution at 
March 4.5 reveals SMC-like extinction in the host galaxy at a 
level of Ay < 0.10, which is significantly lower than for the 
strongly star-forming LBGs. Hence, the extreme faintness of 
the host galaxy indicates a low overall star-formation rate in the 
host galaxy, raising the possibility that the host may be a chem- 
ically less evolved, relatively low-luminosity galaxy containing 
SMC-type dust. We argue that there may be a connection be- 
tween the host galaxy of GRB 00030 1C and DLAs, suggest- 
ing that substantial star-forming activity at high redshift takes 
place in relatively faint galaxies. Future studies of high redshift 
GRBs will further help explore this connection. 

Acknowledgments 

We wish to thank our anonymous referee for providing many 
helpful comments and suggestions. Support for this programme 
by the director of the Nordic Optical Telescope, professor 
Piirola, is much appreciated. We also acknowledge the assis- 
tance given by the ESO service observing team. J. Gorosabel 
acknowledges the receipt of a Marie Curie Research Grant 
from the European Commission. J. Hjorth acknowledges sup- 
port from the Danish Natural Science Research Council (SNF). 
B. Thomsen acknowledges support from the Danish Natural 
Science Research Council for funding the Danish Centre for 
Astrophysics with the HST. G. Bjornsson acknowledges sup- 
port from the Icelandic Research Council and the University 
of Iceland Research Fund. I. Burud is supported by Pole 
d' Attraction Interuniversitaire, P4/05 (SSTC, Belgium). K. 
Hurley acknowledges support for Ulysses operations under 
JPL Contract 958056, for IPN operations under NASA LTSA 
grant NAG5-3500, and for NEAR operations under the NEAR 
Participating Scientist program. We are grateful to R. Gold 
and R. McNutt for their assistance with the NEAR space- 
craft. R. Starr is supported by NASA grant NCC5-380. We 
are indebted to T. Sheets for her excellent work on NEAR 
data reduction. Special thanks also go to the NEAR project 
office for its support of post-launch changes to XGRS soft- 
ware that made these measurements possible. In particular, 
we are grateful to John R. Hayes and Susan E. Schneider 
for writing the GRB software for the XGRS instrument and 
to Stanley B. Cooper and David S. Tillman for making it 
possible to get accurate universal time for the NEAR GRB 
detections. The data presented here have been taken using 
ALFOSC, which is owned by the Instituto de Astrofisica 
de Andalucia (IAA) and operated at the Nordic Optical 
Telescope under agreement between IAA and the NBIfAFG 
of the Astronomical Observatory of Copenhagen. Additionally, 
the availability of the GRB Coordinates Network (GCN) 
and BACODINE services, maintained by Scott Barthelmy, 
is greatly acknowledged. We acknowledge the availabil- 
ity of POSS-II exposures, used in this work; The Second 
Palomar Observatory Sky Survey (POSS-II) was made by 
the California Institute of Technology with funds from the 
National Science Foundation, the National Aeronautics and 



B.L. Jensen et al.: The Afterglow of GRB 000301C 



13 



Space Administration, the National Geographic Society, the 
Sloan Foundation, the Samuel Oschin Foundation, and the 
Eastman Kodak Corporation. 

References 

Allen C. W., "Allen's Astrophysical Quantities", 4th edition 2000, 
Eds. A. N. Cox 

Andersen ML, Hjorth J., Pedersen H. et al., 2000, A&A 364L, 54A 
Berger E. & Frail D. A, 2000, GCN 590 
Berger E., Sari R., Frail D.A., et al., 2000, ApJ 545, 56B 
Bernabei S., Marinoni S., Bartolini C, et al., 2000, GCN 571 
Bertoldi F, 2000, GCN 580 

Blain A.W. & Natarajan P., 2000, MNRAS 312, L39 
Bloom J.S., Diercks A., Djorgovski S.G., et al., 2000, GCN 661 
Bloom J.S., Kulkarni S.R., Djorgovski S.G., et al., 1999, Nature 401, 
453 

Bridger A., Wright G.S., Economou F, et al., 2000, in Ed. H. Lewis, 
Advanced Telescope and Instrumentation Control Software, Proc 
SPIE Vol 4009, in press 
Castro S.M., Diercks A., Djorgovski S.G., et al., 2000, GCN 605 
Castro-Tirado A. & Gorosabel J., 1999, A&AS 138, 449 
Castro-Tirado A., Sokolov V. V., Gorosabel J., et al., 2001, A&A, in 
Press. 

Costa E., et al., 1997, Nature 387, 783 
Dai C.G., Lu T, 1999, ApJ 519, L155 
Feng M., Wang L., Wheeler J. C, 2000, GCN 607 
Fruchter A., Smette A., Gull T, et al., 2000, GCN 627 
Fukugita M., Shimasaku K., Ichikawa T, 1995, PASP 107, 945 
Fynbo J.U., M0ller P., Warren S.J., 1999, MNRAS 305, 849 
Fynbo J.P.U., Jensen B.L., Hjorth J., Pedersen H., Gorosabel J., 2000a, 
GCN 570 

Fynbo J.P.U., Jensen B.L., Hjorth J., Pedersen H, Gorosabel J., 2000b, 
GCN 576 

Fyn bo J.P.U., Jensen B .L., Gorosabel J., et al., 2001a, A&A, In press 

( |astro-ph/0101425| ) 
Fyn bo J.P.U., Goros abel J., Dall T.H., et al., 2001b, A&A, submitted 



(astro-ph/0102158) 



Galama T.J., Vreeswijk P.M., Paradijs J. van, et al., 1998, Nature 395, 
670 

Galama T.J., Tanvir N., Vreeswijk P.M., et al., 1999, ApJ 536, 185 
Garnavich P., Barmby P., Jha S., Stanek K., 2000, GCN 573 
Garnavich P., Loeb_^ J _StaneJ^J£ ; Z ;i _^000 i ^r^ 

ph/0008049) 
Granot J., Piran T, Sari R„ 1999a, ApJ, 513, 679 
Granot J., Piran T, Sari R., 1999b, ApJ, 527, 236 
Haehnelt M.G., Steinmetz M., Rauch M., 2000, ApJ 534, 594 
Harrison FA., Bloom J.S., Frail D.A., et al„ 1999, ApJ 523, L121 
Henden A., 2000, GCN 583 
Henden A., 2001, GCN 961 

Holl and S., Bjornsson G ., Hjorth J., Thomsen B., 2000, A&A, In press 

( astro- ph/0010196| ) 
Hurley K., et al. in Gamma-Ray Bursts, Huntsville 1991, AIP Conf. 

Proc. (AIP New York) 265, p. 3, 1992a 
Hurley K., Sommer M., Atteia J.-L., et al., 1992b, A&AS 92, 401 
Jaunsen A. O., Hjorth J., Bjornsson G, et al., 2001, ApJ 546, 127J 
Klose S., Stecklum B., Fischer O., et al., 2000, GCN 645 
Kobayashi N., Goto M., Hiroshi T, et al., 2000, GCN 577 
Kouveliotou C, et al., 1995, in Proc. 3rd Huntsville Symp. on 

Gamma-Ray Bursts, Ed. C. Kouveliotou, M. F. Briggs & 

G. Fishman, AIP Conf. Proc, New York, p. 42 
Kulkarni S.R., Berger E., Bloom J.S., et al., 2000a, To appear i n Proc. 

of the 5th Huntsville GRB Symposium (|astro-ph/0002168|) 



Kulkarni S.R., Djorgovski S.G., Odewahn S.C., et al., 1999, Nature 
398, 389 

Kulkarni VP, Hill J.M., Schneider G, et al., 2000b, ApJ 536, 36 
MacFadyen A.I., Woosley S.E, 1999, ApJ 524, 262 
Mao S„ Mo H.J., 1999, A&A 339, LI 

Masetti N., Bartolini C, Bernabei S., et al, 2000, A&A 359, L23 
Meszaros P., 1999, to appear in "Black Holes and Gravitational Waves 

- New Eye s in the 21st Century" (eds. T. Nakamura et al.) (a stro- 

ph/99 12546) 

Meszaros P., Rees M.J., 1999, MNRAS 306, 3, L39 
Metzger M.R., Djorgovski S.G., Kulkarni S.R., et al, 1997, Nature 
387, 878 

M0ller P., 2000, The Messenger 99, 3 1 
M0ller P., Warren S.J., 1993, A&A 270, 43 
M0ller P., Warren S.J., 1998, MNRAS 299, 661 
Mukherjee S., Feigelson E.D., Babu G.J., et al., 1998, ApJ 508, 314 
Odewahn S.C., Djorgovski S.G., Kulkarni S.R., et al., 1998, ApJ 509, 
L5 

Paciesas W.S., Meegan C.A., Pendleton G.N., et al., 1999, ApJS 122, 
465 

Panaitescu A., Meszaros P., Rees M.J., 1998, ApJ 503, 314 

van Paradijs L, Groot PL, Galama T, et al., 1997, Nature 386, 686 

Pei Y.C., 1992, ApJ 395, 130 

Piran T, 2000, "Ga mma Ray Bursts: The First Three Minutes" (Ed. 

Juri Poutanen) ( |astro-ph/9907392| ) 
Reichart D.E., 1999, ApJ 521, LI 11 
Rhoads J.E. 1999, ApJ 525, 737 
Rhoads J.E. & Fruchter A.S., 2001, ApJ 546, 1 17R 
Sagar R., Mohan V, Pandey A.K., et al., 2000, Bull. Astron. Society 

of India 28,499 
Sari R., Piran T, Narayan R., 1998, ApJ 497, L17 
Schlegel D.J., Finkbeiner D.P., Davis M., 1998, ApJ 500, 525 
Smette A., Fruchter A., Gull T.R., et al., 2000, GCN 603 
Smette A., Fruc hter A., Gull T.R., et al„ 2001, ApJ, submitted (a stro- 

ph/0007202) 
Smith D.A., Hurley K, Cline T, 2000, GCN 568 
Stanek K.Z., Garnavich P.M., Kaluzny L, Pych W., Thomson I., 1999, 

ApJ 522, L39 

Stecklum B., Klose S., Fischer O., et al., 2000, GCN 572 

Steidel C.C., Adelberger K.L., Dickinson M., Pettini M., 1999, In 

the proceedings of the Xth Rencontres de Blois, "The Birth of 

Galaxies", in press 
Stetson P., 1987, PASP 99, 191S 
Stetson P., 1997, "User's Manual for DAOPHOT II" 
Storrie-Lombardi L.L, Irwin M.J., and McMahon R.G., 1997, 

MNRAS 283, L79 
Totani T, 1997, ApJ 486, L71 
Veillet C, 2000, GCN 575 

Wijers R.A.M.L, Bloom J.S., Bagla J.S., Natarajan P., 1998, MNRAS 
294, LI 3 

Wolfe A.M., Lanzetta K.M., Foltz C.B., Chaffee F.H., 1995, ApJ 454, 
698 

Wolfe A.M., Turnshek D.A., Smith H.E., Cohen R.D., 1986, ApJS 61, 
249 



